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Exact measurement of the second-order correlation function g^^\t) of a light source is essential 
when investigating the photon statistics and the light generation process of the source. For a 
stationary single-mode light source, Mandel Q factor is directly related to For a large 

mean photon number in the mode, the deviation of (0) from unity is so small that even a tiny 
error in measuring would result in an inaccurate Mandel Q. In this work, we have found 

that detector dead time can induce a serious error in g^^\0) and thus in Mandel Q in those cases 
even in a two-detector configuration. Our finding contradicts the conventional understanding that 
detector dead time would not affect g^^\t) in two-detector configurations. Utilizing the cavity-QED 
microlaser, a well-established sub-Poissonian light source, we measured with two different 

types of photodetectors with different dead time. We also introduced prolonged dead time by 
intentionally deleting the photodetection events following a preceding one within a specified time 
interval. We found that the observed Q of the cavity-QED microlaser was underestimated by 19% 
with respect to the dead-time-free Q when its mean photon number was about 600. We derived an 
analytic formula which well explains the behavior of the g^^/O) as a function of the dead time. 

PACS numbers: 42.50.Ar, 42.50.Pq, 42.79.Hp 


I. INTRODUCTION 

The second-order correlation (SOC) function of 

radiation is a key quantity characterizing photon statis¬ 
tics as well as elucidating the underlying light generation 
mechanism. This correlation function is often interpreted 
as being proportional to the probability of measuring a 
photon at time zero and then measuring another photon 
at time t. Since the first observation of SOC of light from 
a star by Hanbury Brown and Twiss (HBT) [T], the mea¬ 
surement techniques for SOC have progressed much as 
high-efRciency photo detectors and fast electronics have 
been developed. 

Significant improvements were made by employing a 
time-to-digital converter or a time digitizer, which pro¬ 
vides a digital representation of the time intervals be¬ 
tween a start photodetection event at one detector and 
multiple stop events at the other detector. From these 
time intervals a histogram of the time delay between 
the start and stop events is obtained. This method is 
called the single-start multi-stop time-to-digital conver- 
sion(SMTDC) In later experiments, a more effi¬ 

cient method which uses all of the arrival time records 
in both start and stop detectors was developed. This 
method is called the multi-start multi-stop time-to- 
digital conversion (MMTDC) compared to the conven¬ 
tional SMTDC. In MMTDC, all of the arrival times at 
both detectors are recorded for a time window Tq much 
longer than the correlation time of a radiation source. 
By software or by using a hardware correlator, we then 
obtain the correlation of all detected photon pairs or 
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more specifically a histogram of time intervals between 
all possible detected photon pairs. The number of detec¬ 
tion events during Tq on a start detector in MMTDC is 
given by Nq = where 77 is the quantum efficiency of 

the detector and is the incident photon flux. Therefore, 
MMTDC is more efficient than SMTDC, which uses only 
one start photon event, by this factor of A^o ^ 1- Owing 
to this high efficiency, MMTDC has been successfully em¬ 
ployed in the first observation of non-classical radiation 
[3 and quantum frequency pulling [5] in the cavity-QED 
microlaser and the spectrum of a single atom localized in 
an optical lattice [7]. 

For accurate measurement of SOC, the effects of detec¬ 
tor characteristics such as detection efficiency and dead 
time have also been investigated, where the latter is a 
time period in which a photo detector becomes blind af¬ 
ter photodetection. Although the SOC function 
of light is independent of detector efficiency, it is appar¬ 
ent that in a single detector configuration detector dead 
time Td can seriously affect the measurement of g^'^'ift). 
Because of the detector dead time, g^‘^\t) is significantly 
reduced for |t| < near the origin. As a result, we 
lose the information on ( 7 O)( 0 ), an important parameter 
directly related to the photon statistics of a stationary 
radiation source as discussed below. No effective way to 
recover the lost information completely has been found 
for a single detector configuration despite many studies 
on this issue [8119]. 

It is often argued that the dead-time deficiency may 
be completely removed in a two-detector configuration 
for stationary light sources. This is based on a simple 
reasoning that two successive photons within the detec¬ 
tor dead time can be resolved if those two photons are 
detected on separate detectors: the first photon is de¬ 
tected on a start detector and the second photon is on a 
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separate stop detector. For non-stationary light sources, 
on the other hand, detector dead time has been shown 
to distort SOC functions as discussed by Choi et al. [10] 
even in a two-detector configuration. 

In this paper, we show that the detector-dead-time 
deficiency cannot be completely removed in the two- 
detector configuration even for stationary light sources. 
In this case, dead time effect comes in as reduction in 
the detected flux due to the missed photons during Td 
compared to the mean waiting time r^, - the mean time 
interval between successive photodetection events. As a 
result, non-negligible distortion occurs in the SOC func¬ 
tion for |t| < Td- The distortion deepens as the 

incident photon flux $ increases and consequently the 
mean waiting time = l/(? 7 $) is reduced to approach 
Td- Such distortion is critical especially for a nonclassi- 
cal light source with a large internal mean photon num¬ 
ber (n). For a stationary single mode light, the relation 
gO)(o) = 1-1- Q/{n) jH] holds with the Mandel Q pa¬ 
rameter bounded between -1 and 0 for nonclassical light. 
For a large mean photon number (n) ^ 1, we then have 
|1 — (0)1 <C 1, so even small distortion by the detec¬ 

tor dead time can cause a large error in determining Q. 
In this work, we first derive a formula quantifying the 
distortion in induced by the dead time and then 

verify its validity in actual experiments with the cavity- 
QED microlaser [^ generating a stationary nonclassical 
radiation. We show that by using the formula we can 
recover (0) and thus the Mandel Q unaffected by the 
detector dead time. 

This paper is organized as follows. In Sec. jnj we dis¬ 
cuss dead time effect on mean photon flux measurement. 
We then extend our discussion to SOC measurement and 
derive a formula to correct the di stor tion introduced by 


the dead time in g^'^^ (0) in Sec. Ill Our experimental 


setup and simulation methods for checking the validity 
of our formula are discussed in Sec. We present the 
experimental and simulation results consistent with our 
theoretical description in Sec. jV] followed by concluding 
remarks in Sec. ED 


II. DEAD TIME EFFECT ON 
PHOTODETECTION FLUX 

Two-detector configuration eliminates the distortion 
due to the missed successive photons on the same detec¬ 
tor by considering two successive photon counting events 
only on separate start and stop detectors. However, there 
still exists another source of distortion coming from the 
reduction in counted photon flux due to the dead time. 


A. light with Poisson photon statistics 

Let us first consider a waiting-time distribution w{t) 
for a detector with a quantum efficiency t] but without 
dead time. If the photon statistics of light is Poisso- 


nian like coherent light, the waiting-time distribution is 
given by a single exponential function: w{t) = 
with (f) = Tyd), the dead-time-free photodetection flux for 
the incident light. The mean waiting time is then 
given by Tu, = In the presence of detector dead 

time, the waiting-time distribution is modified in such a 
way that it vanishes for 0 < t < Td with the rest still 
the same exponential as w(t). When normalized, the 
modified waiting-time distribution w'(t) is nothing but 
w{t — Td). It is then straightforward to show that the 
new mean waiting time t^ is given by 

T'^=Td+Tnj ( 1 ) 

In the presence of the detector dead time, the photode¬ 
tection flux (j)' = 1 /t!^ for the incident light appears to 
be less than the dead-time-free photodetection flux (p by 
the following relation. 


Tw+Td 1 + (pTd 


( 2 ) 


This formula is already derived in the previous works m 
ESj. It has been used to investigate the dead time effect 
on the intensity statistics of scattered light field measured 
with a finite collection aperture HI HD- We can then 
interpret the quantity 4>'/(p as “capture probability” and 
1 — {p’/(p) as “miss probability” in photodetection due to 
the dead time, respectively. 


B. light with non-Poissonian statistics 


If the light source exhibits sub- or super-Poisson pho¬ 
ton statistics, the waiting time distribution is not given 
by a simple exponential function and thus Eq. ([^ is 
no longer valid in general. Eor instance, let us con¬ 
sider light exhibiting sub-Poisson photon statistics with 
its SOC function given by g^^\t) = 1 — with Tc 

the correlation time. Ververk and Orrit [T6| showed that 
the waiting-time distribution is approximately double¬ 
exponential, given by wo(t) ~ for an 

ideal detector of 77 = 1. The relation between p' and p 
would then be quite different from Eq. ([^. 

In general, the capture probability can be written in 
terms of the detector dead time as 


where x = pTd and the coefficient is given by 


an 


1 d^{P/P’) 


n\ dx'^ 


x—0 


(4) 


depending on the specific waiting time distribution of the 
system under consideration. For the above waiting time 
distribution WQ{t), the lowest non-vanishing coefficient is 
02 = [2(1 — pTc)pTc]~^ under the condition Td < Tc- 
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If the mean internal photon number (n) of a source is 
much larger than \Q\ regardless of its photon statistics, 
the SOC function is close to that of coherent light, i.e., 
^ 9,nd the corresponding waiting-time dis¬ 
tribution is approximately single-exponential. An exam¬ 
ple is the cavity-QED microlaser, where (n) ^ 10^ — 10^ 
and — 1 < Q < 1. We can then use Eq. for considera¬ 
tion of a detector dead time on the photodetection flux. 
The precise condition for the validity of this approxima¬ 
tion is |1 — ^ 1 as shown in Appendix. 

Ververk and Orrit neglected the higher-order correla¬ 
tion effects in deriving the above waiting-time distribu¬ 
tion. More general discussion including the higher-order 
correlations is presented in Ref. El US]. We also neglect 
the higher-order correlations (0,..., 0, t) with n > 3 
in this paper. It is because that the higher-order corre¬ 
lations are not much different from SOC, unlike thermal 
light sources m, in the cavity-QED microlaser, where 
the condition |1 — <C 1 is well satisfied. 


III. DEAD TIME EFFECT ON SECOND-ORDER 
CORRELATION MEASUREMENT 


In a two-detector configuration, the intensity correla¬ 
tion {Ist{t)Isp{t + t')) is understood as a joint probability 
of photodetection at time t on a start detector and at 
time t + t' on a stop detector. The intensity operator for 
the start(stop) port is denoted as Ist{Isp)- For station¬ 
ary light, the correlation does not depend on t and thus 
it can be replaced with a hxed time. If we define JV{t) to 
be the actual number of photon pairs comprised of one 
photon incident on the start detector and another on the 
stop detector with a time delay t, the normalized SOC 
function (t) in this conhguration can be expressed as 

(2)/i'| _ {ISt jo)Ispjt)) _ A/'(t) _ N{t) , , 

^ ^ ^ " (Istmilspit)) ~ AA(oo) " iV(oo) ’ ^ ^ 


start I 

stop I 


y 
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FIG. 1. Illustration of photodetection events in a HBT-type 
two-detector configuration. The number of photon pairs with 
a zero time delay (circled on the left) and that with a very 
long time delay (circled on the right) are needed in evaluating 


probable by the capture probability 4)'/(j) = (1 + (j^Td) 
so the denominator N{oo) should be replaced by 


A/'(oo) —>• 


A/'(oo) 

(1 -h (j)stTd) (1 + 4>spTd) ’ 


( 6 ) 


where 4‘st{4‘sp) is the dead-time-free photodetection flux 
on the start (stop) detector. 

For zero time delay, on the other hand, the photo flux 
on each detector is further modified by a factor g^^\0). 
It is because the probability of having a photon on one 
detector with another photon on the other detector at 
the same time is proportional to g^‘^\0). Including this 
effect, the photo detection flux (j) is replaced by g^‘^\0)(j) 
for each detector. The numerator A/’(0) should then be 
replaced by 


AA(0) ^ 


AA(0) 


{1 + g^^H0)(l)stTd) + g^'^K0)(4spTd) 
As a result, the observed SOC will be 


(7) 


9-gl(0) = 9(^>(0) ‘ + (’•‘0 ‘ . (8) 

1 -I- [OjCpstTd 1 + g^^> [0)(j)spTd 


When the photodetection flux on each detector is as 
low as (jjTd 1, Eq. Q is further approximated as 


where N{t) is the number of photodetection pairs with 
a time delay t in the absence of detector dead time, so 
N{t) = Equation ([^ shows that SOC does not 

depend on the detector efficiency. 

From now on, let us concentrate on g(^)(0), a parame¬ 
ter directly related to the Mandel Q of the internal field 
of the source. In evaluating Eq. (^, the 

numerator A/'(0) is obtained by counting events like the 
circled one on the left in Fig. whereas the denominator 
N{oo) is obtained from the events like the one circled 
on the right. Our interest is then how the detector dead 
time affects such counting events. We restrict ourselves 
to the case of \1 — g^‘^\t)\ 1, i.e., to the case where the 

waiting-time distribution is near single-exponential and 
thus we can still use Eq. ([^. 

Let us first consider the time delay t much larger than 
a correlation time. In this case, we can neglect the corre¬ 
lation between photons in each pair. Because of the de¬ 
tector dead time, each photo counting event is then less 


(0) ~ 5^^) (0) 11 -b (0) {4>st + (4sp)Td } 

or in terms of Mandel Q 

Q' 1 - (0)((/)st -b (l>sp)Td 


(9) 


( 10 ) 

This approximation shows that the dead time effect be¬ 
comes important as the photodetection flux increases for 
a given dead time. 

For an arbitrary waiting time distribution, the capture 
probability of photodetection under detector dead time 
can be written as that in Eq. ([^. If the lowest non¬ 
vanishing term in the denominator is of the n-th order, 
Eq. (10) is then replaced with 

} ■ (H) 

with an is given by Eq. Q. 
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IV. EXPERIMENT AND DEAD-TIME 
SIMULATION 


A. Counter Electronics 

The counter electronics used in the present experi¬ 
ment is an improved version of the one used by Choi 
et al. [smo]. In the former system, two counter/timing 
boards were separately installed in two personal com¬ 
puters (PC’s) in order to avoid interchannel crosstalk. 
Another PC was used to trigger those counter/timing 
boards and to control overall measurement sequence. In 
this configuration, each counter/timing board, once trig¬ 
gered, records photodetection times based on its own in¬ 
ternal clock, and thus the clocks in those counting chan¬ 
nels were not synchronized. 

Clock synchronization is realized in the present setup 
by employing a counter board (National Instruments NI- 
7813R) equipped with a held programmable gate array 
(FPGA). By programming the FPGA we have imple¬ 
mented multiple counting channels without crosstalk in a 
single board. Those counting channels are perfectly syn¬ 
chronized at a clock speed of 125 MHz while internally 
triggered. The resulting time resolution is 8 ns, improved 
from 12.5 ns of the former system. 

Moreover, with the present setup using an FPGA we 
have eliminated the counter-board-related dead-time ef¬ 
fect observed in the previous setup. When a photode¬ 
tection event occurs, the counter board in the previous 
setup saves the event time measured in clock period in 
an onboard register hrst and then it is transferred to the 
computer memory through direct memory access (DMA). 
A problem arises when the next photodetection event oc¬ 
curs before this transfer is completed: the new event is 
simply ignored. As a result, the SOC in a single channel 
configuration exhibits a partial dead-time effect as shown 
in Fig.[2];a), where the SOG function obtained for coher¬ 
ent light is plotted. The narrow perfect dip in the range 
of |<| < 50 ns is due to the detector dead whereas the wide 
partial dips in the range of 100 ns < |t| < 250 ns arises 
from the above data loss at the counter board. We call 
this the counter dead-time effect. In the present setup us¬ 
ing an FPGA, the onboard hrst-in-first-out (FIFO) mem¬ 
ory acts as a large buffer in the DMA transfer and thus 
can eliminate the above counter dead-time effect. Its per¬ 
formance in a single-channel SOG function measurement 
is shown in Fig. |^b), where only the detector dead time 
effect is noticed in the range of |t| < 20 ns from without 
any partial dips due to the counter dead time. Glear iso¬ 
lation of the detector dead time effect as in Fig. Hb) in 
fact enables us to correct the SOG function against the 
detector dead time in Sec. El 



time delay (gs) 


FIG. 2. Counter-board-related dead-time effect in the ob¬ 
served in a single detector configuration, (a) With the 

previous counter/timing board with an onboard register, (b) 
With the new FPGA counter board with an onboard FIFO 
memory. Poissonian background light was used as a light 
source. The partial dips in the range of 100 ns < |t| < 250 ns 
are due to the data transfer loss between the counter/timing 
board and a control computer. 


dead time of 50 ns (Perkin Elmer SPGM-AQR-12) and 
it will be referred as SPGM-S (slow). The other has a 
shorter dead time of 21 ns (Excelitas SPGM-AQRH-12) 
and it will be referred as SPGM-E (fast). They have the 
same characteristics except for the dead time and out¬ 
put voltage specification. The output voltage pulse of 
SPGM-F is not compatible with our counter board (NI- 
7813R). A homemade pulse stretcher made of fast logic 
gates is used between them: the output pulse width of 7 
ns of SPCM-E is extended to 12.5 ns with an enhanced 
peak height for the counter board as shown in Eig. 


B. Single-Photon Counting Detectors 


C. Experiment 


Two different models of single-photon counting mod¬ 
ules (SPGM’s) are used in our experiment. One has a 


Our experimental schematic is depicted in Fig. The 
basic physical principles and apparatus to generate sub- 
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FIG. 3. The shape of a pulse (black) from a SPCM-F 
and that of a pulse (blue) processed by a homemade pulse 
stretcher. In the pulse stretcher, the original pulse and its 
delayed pulse are added in time with a small gain by using a 
OR gate. A delay of a few nanoseconds is due to the intrinsic 
time delay of the gate chip. 



cavity 

locking 

beam 



magnetic field 


FIG. 4. Experimental schematic. The cavity-QED micro- 
laser is pumped by a beam of barium atoms prepared in the 
excited state [^. The SOG of the output is measured in a two- 
detector configuration. Flippable mirrors (FM’s) are used to 
select a desired SPCM pair, SPCM-F or SPCM-S, while keep 
the other experimental conditions unchanged. 


Poisson light with the cavity-QED microlaser is the same 
as in the previous work by Choi et al. . In order to fa¬ 
cilitate the switching between the two types of detectors 
with different dead times, flippable mirrors are used to 
provide a choice of detectors while preserving the other 
experimental conditions. 

The SPCM manufacturers provide empirical counting 
correction factors in the instruction manual |20) up to 
photodetection flux of 2.5 x 10^ counts/sec or 25 Mcps 
(Mega counts per second). These correction factors start 
to deviate from that in Eq. ([^ at photon flux of 4 Mcps 
by unexplained reasons. In order to avoid this inconsis¬ 
tency at high photon flux, we attenuate the photon flux 
to keep the photodetection flux under 3 Mps, where the 
manufacturers’ correction factors well agree with Eq. ([^. 

We have also measured the detector dead times from 
the actual waiting-time distributions. With a combina¬ 
tion of two SPCM-F’s and the EPGA counting board. 


(a) 

(i) • 


(ii) #“ 


(iii) 


(iv)#“ 


• • 
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FIG. 5. Over-deletion problem and its effect on the waiting¬ 
time distribution. In (a), the red -filled circles represent pho¬ 
todetection events in sequence whereas the pink-filled circles 
indicate missed events due to an intrinsic or a prolonged dead 
time. (i). Dead-time-free case. (ii). With an intrinsic dead 
time, whose length is indicated by a solid arrow. The third 
event is missed due to the intrinsic dead time. (iii). With 
a prolong dead time. The prolonged part from the intrinsic 
dead time is marked by a dashed arrow. The second event 
within the prolonged dead time is deleted, (iv). With an 
intrinsic dead time as long as the prolong dead time in (iii). 
Differently from case (iii), the third event is detected. There¬ 
fore, over-deletion occurs in (iii) compared to the case of a real 
dead time, (b) Waiting time distribution with a 80-ns pro¬ 
longed dead time from a 28-ns intrinsic dead time. A small 
distortion, marked by a blue arrow, occurs around 80 ns, but 
it is so small that the waiting time distribution can still be 
approximated by a single exponential function. 


the dead time extracted from the waiting-time distribu¬ 
tion was 28 ns. For SPCM-S, the observed dead time was 
56 ns. Detection bin time was 8 ns for both cases. 

We neglect the effect of after-pulsing in our measure¬ 
ment because the probability of after-pulsing is only 
0.3% per real photodetection according to the detector 
manual. Furthermore, after-pulses at separate detectors 
are perfectly uncorrelated, and thus they contribute a 
Poissonian background in SOC measurement. Another 
Poissonian background associated with the detector dark 
counts is less than 0.1% of the actual photon flux in our 
experiment, so it does not affect our correlation mea¬ 
surement either. If a Poissonian background count rate 
were comparable to the actual photon flux, detector dead 
time could make a huge difference in measuring SOC in 
a different way as previously studied in Ref. [21]. 


D. Simulating Prolonged Detector Dead Times 

In order to investigate the detector dead time effect 
systematically in experiment, it is desired to have as 
many detectors with the same characteristics but with 
different dead times as possible. In reality, we have a 
limited number of detectors. We have two SPCM-F’s 
with a mean dead time of 28-ns and two SPCM-S’s with 
a 56-ns dead time. 

To overcome this limitation, we have simulated addi¬ 
tional dead times for a given detector by deleting the 


































6 


subsequent photodetection records within an extended 
period beyond the actual dead time after any photode¬ 
tection event. Those extended periods serve as prolonged 
dead times. 

However, the effect of a prolonged dead time is not ex¬ 
actly the same as that of a real dead time with an equal 
magnitude. It is because with a prolong dead time some 
photodetection events are lost which would be detected 
with a real dead time with the same magnitude as illus¬ 
trated in Fig. [^a). This over-deletion of counts leads 
to a distortion in the waiting time distribution function 
near the prolonged dead time as indicated in Fig. [^b). 
Nonetheless, the distortion is not large enough to affect 
the overall shape of the waiting time distribution, indi¬ 
cating the over-deletion rarely occurs. Therefore, we can 
utilize the prolonged dead times for systematic investiga¬ 
tion of the dead time effect in the next section. 


V. RESULTS AND DISCUSSIONS 

We have measured (0) (in the presence of the dead 
time effect) of the cavity-QED microlaser output by us¬ 
ing the setup depicted in Fig. with two different sets of 
detectors. The results are shown in Fig. [^a), where the 
red star indicates the result with SPCM-F’s and the pur¬ 
ple square corresponds to the result with SPCM-S’s. The 
dead-time-corrected photodetection flux on each detector 
was 2.6 Mcps for SPCM-F and 3.3 Mcps for SPCM-S. 
Also plotted in Fig. |6l;a) as black dots are the results 
obtained with prolong dead times as discussed in Sec. 
HvlC on SPCM-F’s of 28-ns dead time. A solid curve is 
a theoretical fit by Eq. Q, well agreeing with the results 
with actual and prolong dead times. The only fitting 
parameter is which appears as a vertical axis off¬ 

set corresponding to zero dead time. The dead-time-free 
fluxes needed in Eq. (|^ are obtained from the observed 
photodetection fluxes by using Eq. (§. 

The smallest g'^^\0) — 1 measured with SPCM-E is 
-710±60 ppm (part per million), corresponding to Man- 
del Q of -0.43±0.04, whereas the dead-time-corrected 
(0) — 1 obtained from the fitting is -850±60 ppm and 
thus the actual Mandel Q of the microlaser output is - 
0.51±0.04. It is noteworthy that the smallest dead time 
effect with our best detector still amounts to a consider¬ 
able distortion (0.51/0.43-1=19%) in the Mandel Q mea¬ 
surement. The dead time correction by Eq. Q is thus 
essential for accurate photon statistics measurement. 

In Fig.§b), we examine the dependence of (0) on 
photodetection flux. The black and red squares repre¬ 
sent the results of ( 7 ^^^( 0 ) measurement with the same 
detectors (SPCM-S) but with a reduce flux for the red 
one. The solid curves are theoretical fits by Eq. Q with 
(0) as a fitting parameter individually. The fit curves, 
although independently done, almost meet at zero dead 
time. The discrepancy is less than the experimental er¬ 
ror. 

The observed dependence g'^^^O) on the photodetec- 



FIG. 6. Dead time effect on the SOC function of the cavity- 
QED microlaser. (a) Dependence of the observed SOC on 
intrinsic and prolonged dead times. Red star is the observed 
with SPCM-F’s of a 28-ns dead time whereas the pur¬ 
ple square is that with SPCM-S’s of a 56-ns dead time. Black 
dots are the results obtained with prolonged dead times ap¬ 
plied to the photodetection records corresponding to the data 
marked by the red star. Blue line is a theoretical fit given 
by Eq. (|^ under the condition that it should go through the 
red dot. (b) Dependence of the observed SOC on photode¬ 
tection flux. Black(red) square is the observed g'^^^(O) with 
SPCM-S’s when the photo fluxes on start and stop detec¬ 
tors are 2.75(1.94) Mcps and 2.23(1.60) Mcps, respectively. 
Black(red) solid curve is a theoretical fit given by Eq. ^ un¬ 
der the condition that it should go through the black(red) 
square. Error bars represent the fitting errors in obtaining 
values from the SOC data. Typical fitting errors for 
simulated data with prolonged dead times are normally 10% 
of |g^^^(0) ~ l|' The error bar for the black square in (b) is 
smaller than than size of the square. 


tion flux may appear contradictory to the general view 
that the SOC function does not depend on random miss 
of incident photons. An example is detector efficiency. 
It should be noted, however, that the missing of incident 
photons due to the detector dead time is not random 
miss at all. The missing occurs only immediately after 
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a successful photodetection event. In other words, the 
missing event is correlated with the success event with 
the detector dead time as a correlation time. So, there is 
no contradiction. 

One way to avoid the dead time effects discussed so far 
is to keep the photodetection flux low enough to make 
the mean waiting time very much larger than the dead 
time, The distortion in SOC will then be 

negligible as shown in Eq. Q. This approach does not 
work, however, when the measurement time Tq for SOC 
is practically limited. For instance, a photon source may 
have a finite operating time and thus Tq is limited. In 
order to resolve the feature in near the origin for 

Mandel Q measurement, the signal-to-noise ratio has to 
be larger than |1 — For a bin time tb «C Tc, we 

then have to satisfy y// (ro/th)|l — > 1 

or Tyj < -s/Tot&ll — This requirement set a upper 

bound for . If this upper bound is not much larger than 
the detector dead time, we cannot avoid the dead time 
effect and thus have to rely on our correction formula. 
In the original HBT experiment, the photon flux from 
a distant star were extremely low. Consequently, the 
waiting time was much larger than the detector dead 
time and thus its effect was negligibly small. 


VI. CONCLUSION 

We have found that detector dead time can induce a 
significant error in the measurement of it) of a sta¬ 
tionary light source even in two detectors configuration 
when the internal mean photon number is much larger 
than unity. This finding contradicts the conventional 
understanding that detector dead time would not affect 
(t) in two-detector configurations. In experiment, we 
employed the cavity-QED microlaser for a sub-Poissonian 
light source and measured (0) with two different types 
of photodetectors with different dead time. The ob¬ 
served Q of the cavity-QED microlaser was underesti¬ 
mated as much as 19% with respect to the dead-time-free 
Q even when we used single-photon counting modules 
with the shortest dead time available. We also simu¬ 


lated prolonged dead times by intentionally deleting the 
photodetection events following a preceding one. The 
observed values of g^‘^'> (0) for various real and prolonged 
dead time were well explained by our analytic formula. 
Dead-time-free thus Mandel Q of a stationary 

light source can be obtained with our correction formula. 

The present work is limited to the case of a large 
internal mean photon number so as to neglect the 
higher-order correlation. By including the higher-order 
correlation, one may obtain a more general formula for 
the detector dead time effect on 


Appendix: Conditions for approximating the waiting 
time distribntion as a single exponential 

Let us consider the Laplace transforms W{s) and G{s) 
of a waiting-time distribution wit) and a SOC function 
g^^^t), respectively. According to Ref. [T5], they are re¬ 
lated as 


W(s)=</.G(s)/(l + </)G(s)). (A.l) 


For the cavity-QED microlaser operating at a high mean 
photon number, gi'^i (t) can be written as 

g(2)(t) = (A.2) 


Then G(s) is readily given by 


G(s) = - + ^, 
s s + a 


(A.3) 


where a = IjTc- Inserting this expression in Eq. (A.l) 


and taking inverse Laplace transform, we obtain the 
waiting-time distribution of the microlaser as 


w{t) = (/)£ ^ 


(1 — /3)s -I- a 

s2 -I- [a -I- (j){l - /3)] s -I- (jia 


(A.4) 


The above waiting-time distribution can be approxi¬ 
mated by a single exponential if /3 ^ 1, which can 

be expressed as |1 — ^ 1- 
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